Quantum entanglement and teleportation in pulsed cavity-optomechanics 
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Entangling a mechanical oscillator with an optical mode is an enticing and yet a very challenging 
goal in cavity-optomechanics. Here we consider a pulsed scheme to create EPR-type entanglement 
between a travelling-wave light pulse and a mechanical oscillator. The entanglement can be verified 
unambiguously by a pump-probe sequence of pulses. In contrast to schemes that work in a steady- 
state regime under continuous wave drive this protocol is not subject to stability requirements that 
normally limit the strength of achievable entanglement. We investigate the protocol's performance 
under realistic conditions, including mechanical decoherence, in full detail. We discuss the relevance 
of a high mechanical Q ■ /-product for entanglement creation and provide a quantitative statement 
on which magnitude of the Q ■ /-product is necessary for a successful realization of the scheme. We 
determine the optimal parameter regime for its operation, and show it to work in current state-of- 
the-art systems. 



I. INTRODUCTION 

In optomechanical systems a cavity mode can be 
strongly coupled to a high-quality mechanical oscilla- 
tor via radiation pressure or dipole gradient forces [1—3]. 
Quantum effects [4-6] are starting to play an increasingly 
important role: In the microwave regime ground-state 
cooling via laser-cooling techniques [7], strong coupling 
[8, 9] and coherent control of single-phonon excitations 
[9] have been successfully achieved. In the optical regime 
cooling to the quantum ground state [10] and effects of 
strong coupling [11] have been demonstrated in recent 
experiments. It is as yet an outstanding goal to ob- 
serve genuine quantum effects such as entanglement [12] 
at macroscopic length and mass scales. 

Entanglement of a mechanical oscillator with light has 
been predicted in a number of theoretical studies [13-23] 
and would be an intriguing demonstration of optomc- 
chanics in the quantum regime. These studies, as well 
as similar ones investigating entanglement among sev- 
eral mechanical oscillators [24-32] , explore entanglement 
in the steady-state regime. In this regime the optome- 
chanical system is driven by one or more continuous- wave 
light fields and settles into a stationary state, for which 
the interplay of optomechanical coupling, cavity decay, 
damping of the mechanical oscillator, and thermal noise 
forces may remarkably give rise to persistent entangle- 
ment between the intracavity field and the mechanical 
oscillator. 

Entanglement in the steady-state regime shows two 
main characteristic features: Firstly, entanglement 
reaches a maximal value when the system is driven close 
to a point of dynamical instability. Vice versa, the lim- 
its on the strength of entanglement achievable in proto- 
cols working in the steady-state regime are set by the 
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very conditions guaranteeing a dynamically stable, sta- 
tionary state. Recent studies indicate that these limi- 
tations can even become rather restrictive when a finite 
laser linewidth is taken into account [20, 22]. Secondly, 
the verification of entanglement between the intracav- 
ity field and the moving mirror has to be performed via 
measurements on the outcoupled light leaving the op- 
tomechanical system. Ultimately only correlations be- 
tween modes of the light field are measured, from which 
any entanglement involving the mechanical oscillator has 
to be inferred. However, due to the curious feature of 
quantum correlations that "no entanglement is necessary 
to distribute entanglement" [33, 34], this sort of infer- 
ence is in general a delicate issue. It is unambiguously 
only possible under additional assumptions regarding the 
particular dynamics (i. e., the system's Hamiltonian) and 
structure of the steady state [13]. 

An alternative approach to achieving optomechanical 
entanglement is to work in the pulsed regime, where en- 
tanglement is created and verified with two subsequent 
pulses of light. This strategy has first been developed 
in the context of atomic ensembles [35] and was recently 
considered for systems employing levitated microspheres 
trapped in an optical cavity [36]. A pulsed scheme does 
not rely on the existence of a stable steady state, which 
provides us with the benefit that entanglement is not lim- 
ited by stability requirements. The temporal ordering of 
the pulses excludes the possibility of distributing entan- 
glement without using entanglement [33, 34] such that it 
provides a direct (i.e., without additional assumptions) 
and unambiguous test of entanglement. Similar protocols 
have also been discussed for micromirrors in free space 
(i. e., without the use of an optical cavity) [37, 38]. 

In this article we provide a complete treatment of a 
protocol for the generation and verification of optome- 
chanical entanglement using pulsed light. In addition to 
an idealized scenario, which was briefly discussed in [36], 
we include in our description the full dynamics of the 
optomechanical system. Our derivation provides an ex- 
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haustive discussion of imperfections, how they affect the 
performance of the protocol and how these effects can 
be minimised. Most prominently, we (perturbatively) in- 
clude thermal decoherence of the mechanical system and 
find that a high Q- /-product (quality factor times the fre- 
quency of the mechanical oscillator) plays a crucial role 
in the creation of optomechanical entanglement. More 
specifically, we find that the relation Q ■ / ^> k^T/h 
(where T is the temperature of the mechanical environ- 
ment) has to hold, which is a general and often very 
stringent condition to observe quantum effects in optome- 
chanical systems [39]. A large Q ■ /-product is thus one 
of the most important aspects to consider in the design 
of novel high-quality mechanical resonators. 

To further explore the effect of imperfections on the 
protocol, we optimise the amount of created entangle- 
ment with respect to key experimental parameters and 
present specific values for two existing optomechanical 
systems. We find that creation of entanglement is possi- 
ble in a parameter regime which is realistic yet challeng- 
ing for current state-of-the-art setups. Very importantly, 
our treatment also provides a quantitative statement on 
what magnitude of the Q ■ /-product is necessary in or- 
der to successfully create entanglement for a given bath 
temperature T. These findings provide a general under- 
standing of the requirements to observe quantum effects 
in optomechanical systems and represent essential infor- 
mation for the material development and for the further 
design of future optomechanical structures. 

We finally note that the quantum state created in 
this protocol exhibits a type of entanglement known as 
Einstein-Podolsky-Rosen (EPR) entanglement [40] be- 
tween the mechanical oscillator and the light pulse [36]. 
It thus provides the canonical resource for quantum in- 
formation protocols involving continuous variable (CV) 
systems [41, 42]. We give a detailed description of how 
optomechanical EPR entanglement can be used for the 
teleportation of the state of a propagating light pulse 
onto a mechanical oscillator. 

Note that other promising perspectives for pulsed op- 
tomechanics have recently been discussed also in [43- 
47] , albeit in a very different parameter regime employ- 
ing light pulses which are short on the time scale of a 
mechanical oscillation. The importance of temporal or- 
dering in the verification of optomechanical entanglement 
has also been pointed out in [15]. 

The paper is organized as follows: Section II A contains 
the main results of this work. We first describe entangle- 
ment creation under idealized circumstances and outline 
a way to verify it unambiguously. Additionally, we show 
how it can be used as a resource for CV teleportation. 
In Sec. II B we analyze the influence of imperfections on 
the protocol's performance and find the optimal param- 
eter regime for maximal entanglement. Section III gives 
a detailed description of the full system dynamics. The 
Appendix contains a short derivation of the effective sys- 
tem Hamiltonian in the pulsed regime. 
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FIG. 1. (Color online) Schematic of the system and the tele- 
portation protocol: (a) A blue detuned light pulse (A) is en- 
tangled with the mirror (B). (b) A second light pulse (V) is 
prepared in the input state and interferes with A on a beam 
splitter. Two homodyne detectors measure P° ut + X v and 
X° +P V , yielding outcomes mx and mp, respectively. Feed- 
back is applied by displacing the mirror state in phase space 
by a unitary transformation flxjmx) DpJyUip). (c) To ver- 
ify the success of the protocol, the mirror state is coherently 
transferred to a red detuned laser pulse and a generalised 
quadrature X[{0) = X{° ut cos 9 + jy° ut sinf? is measured. Re- 
peating steps (a)-(c) for the same input state but for different 
phases 8 yields a reconstruction of the mirror's quantum state. 

II. CENTRAL RESULTS 
A. Motivation for the pulsed scheme 

1. Cavity optomechanical system 

Let us consider an optomechanical cavity in a Fabry- 
Perot-type setup (Fig. 1), with mechanical oscillation fre- 
quency uj m , mechanical dissipation rate 7, optical reso- 
nance frequency uj c and cavity decay rate k. A light pulse 
of duration r and carrier frequency oj\ impinges on the 
cavity and interacts with the oscillatory mirror mode via 
radiation pressure. 

In a frame rotating with the laser frequency, the system 
is described by the (effective) Hamiltonian [48] 

H = w m <4 a ro + A c aja c + g (a m + a*,) (a c + oj) (1) 

where a m and a c are annihilation operators of the me- 
chanical and optical mode, respectively. The condi- 
tions under which (1) is valid are discussed in detail 
in Sec. Ill A. The first two terms give the energy of 
the mechanical oscillator and the cavity field, where 
A c = uj c — u)\ is the detuning of the laser drive with respect 
to the cavity resonance. The last term describes the lin- 
earised optomechanical coupling (with coupling constant 
g) via radiation pressure x m x c cx (a m + a^) (a c + aj) = 
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(a m a| + aj n a c ) + (a m a c + a m al) , which can be decom- 
posed into two terms [4, 25, 37, 38]: a beam-splitter-like 
interaction (the first term) and a two-mode-squeezing in- 
teraction (the second term). The former can be used to 
cool the mirror as well as to generate a state swap be- 
tween the mechanical and the optical mode, while the 
latter term describes the optomechanical analogue to the 
optical down-conversion process in an optical paramet- 
ric amplifier and is known to create entanglement from 
coherent input states [4, 49]. 

We assume the pulse to approximately be a flat-top 
pulse, which has a constant amplitude for the largest 
part, but possesses a smooth head and tail (see the Ap- 
pendix). The coupling constant g is then given by 
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with ./Vph being the number of photons in the pulse (see 
Sec. Ill A) . The single-photon coupling constant g is de- 
fined by go = uj c xq/L, where xq is the size of the zero- 
point-motion of the mechanical oscillator and L is the 
cavity length. It is possible to make a single one of the 
interaction terms dominant by tuning the laser such, that 
one of its motional sidebands L0\±uj m is resonant with the 
cavity, where for blue detuning the resonant scattering to 
the lower (Stokes) sideband (u) c — uj\— uj m ) enhances the 
down-conversion interaction, while for red detuning the 
resonant scattering to the upper (anti-Stokes) sideband 
(cj c = uj\ + LJ m ) enhances the beam-splitter interaction 
[4]. In this pump-probe scheme we make use of both 
dynamics separately: Pulses tuned to the blue side of 
the cavity resonance are applied to create entanglement, 
while pulses on the red side are later used to read out 
the final mirror state. A similar separation of Stokes and 
anti-Stokes sideband was suggested in [37, 38] by select- 
ing different angles of reflection of a light pulse scattered 
from a vibrating mirror in free space. 

The full system dynamics, including the dissipative 
coupling of the mirror and the cavity decay, are de- 
scribed by quantum Langevin equations [50], which de- 
termine the time evolution of the corresponding opera- 
tors x m = (a m + <4j/\/2, Pm = ~K a m - «m)/%/2 and 
a c , aj. They read 

X m = ^mPm, (3a) 

p m = -u m x m - 7p m - V2g (a c + aj) - \/2j f, (3b) 
a c = — (iA c + n)a c — iV2gx m — V2/ta in , (3c) 

where we introduced the (self-adjoint) Brownian stochas- 
tic force /, and quantum noise a in entering the cav- 
ity from the electromagnetic environment. Both a in 
and — in the high-temperature limit — / are assumed to 
be Markovian. Their correlation functions are thus given 
by (a in (t)a\ n (t')) = 5(t — t') (in the optical vacuum state) 
and (f(t)f(f) + f(f)f(t)) = (2n + l)S(t - t') (in a ther- 
mal state of the mechanics) [50] . 



2. Creation of optomechanical entanglement 

In this section we impose the following conditions on 
the system's parameters. Firstly, we drive the cavity with 
a blue detuned laser pulse (A c = — uj m ) and assume to 
work in the resolved-sideband regime (k <C u> m ) to en- 
hance the down-conversion dynamics. Note that in this 
regime a stable steady state only exists for very weak op- 
tomechanical coupling [51], which poses a fundamental 
limit to the amount of entanglement that can be created 
in a continuous-wave scheme [14]. In contrast, a pulsed 
scheme does not suffer from these instability issues. In 
fact, it is easy to check by integrating the full dynamics 
(see Sec. Ill B) up to time r, that working in this particu- 
lar regime yields maximal entanglement, which increases 
with increasing sideband resolution u m / k. Secondly we 
assume a weak optomechanical coupling j«k, such that 
only first-order interactions of photons with the mechan- 
ics contribute. This minimises pulse distortion and sim- 
plifies the experimental realization of the protocol. Taken 
together, the conditions g -C k <C w m allow us to invoke 
the rotating- wave approximation (RWA), which amounts 
to neglecting the beam-splitter term in (1). Also, we ne- 
glect mechanical decoherence effects in this section. We 
emphasise that this approximation is justified as long as 
the total duration of the protocol is short compared to 
the effective mechanical decoherence time 1/771, where 7 
is the mechanical damping rate and n is the thermal oc- 
cupation of the corresponding bath. Corrections to this 
simplified model — including the treatment of mechani- 
cal decoherence and dynamics beyond the RWA — will be 
addressed in Sec. II B. 

Based on the assumptions above we can now simplify 
equations (3). For convenience we go into a frame rotat- 
ing with w m by substituting a c — > a c e IWm *, 



and a m — > a m e _I " m . Note that in this picture the cen- 
tral frequency of a in is located at ui\ — uj m = w c . In the 
RWA the Langevin equations then simplify to 
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In the limit g k we can use an adiabatic solution for 
the cavity mode and we therefore find 
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a m {t) a e G *a m (0) +iV2Ge Gt / ds e" Gs aUs), (5b) 



where we defined G — g 2 / k. Equation (5b) shows that 
the mirror motion gets correlated to a light mode of cen- 
tral frequency uj\ — uj m (which coincides with the cavity 
resonance frequency lu c ) with an exponentially shaped 
envelope a- m (t) cx e~ Gt . Using the standard cavity input- 
output relations a out = a in + \/2k a c allows us to define 
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a set of normalised temporal light modes 

A„ = y 1 _ c _ 2GT J dt e" Gt a in (t), (6a) 

2G /" r 
e2GT _ 1 J dt e Gt a out (i), (6b) 

which obey the canonical commutation relations 
[AjjAj] = 1. Together with the definitions B m = a m (0) 
and .Bout = a m ( T ) we & rr i ve &t the following expressions, 
which relate the mechanical and optical mode at the end 
of the pulse t = r 

A out = - e GT A„ - We 2GT -lBl, (7a) 

Bout = e GT B in + iVe 2G - -l4 n . (7b) 

By expressing equations (7) in terms of quadratures 
X % m = (Bi + B\)/y/2 and X( = {A, + Al)/\/2, where 
i G {in, out}, and their corresponding conjugate vari- 
ables, we can calculate the so-called EPR variance Aepr 
of the state after the interaction. For light initially in 
vacuum (AX™) 2 = (AP" 1 ) 2 = \ and the mirror in a 
thermal state (AX™) 2 = (AB™) 2 = n + 1, the state is 
entangled iff [52] 

A EPR = [A(XT + B° ut )] 2 + [A(B° ut + X? ut )] 2 

( , \ 2 ( g ) 

= 2(n + l) (e''-v / e 2 ' -lj < 2, 

where r = Gr is the squeezing parameter and no the ini- 
tial occupation number of the mechanical oscillator. Note 
that in the limit of large squeezing r>lwe find that the 
variance Aepr w (n + 1) e~ 2r /2 is suppressed exponen- 
tially, which shows that the created state asymptotically 
approximates an EPR state. Therefore, this state can be 
readily used to conduct optomechanical teleportation as 
described in Sec. II A 4. 

Rearranging (8) , we find that the state is entangled as 
long as 

1 / («o+2)^\ n ^oo 1 , . 

r > rn = - In — r ~ -lnnn. ( 9 ) 

This illustrates that in our scheme the requirement on 
the strength of the effective optomechanical interaction, 

2 

as quantified by the parameter r = scales logarith- 
mically with the initial occupation number no of the me- 
chanical oscillator. This tremendously eases the proto- 
col's experimental realization, as neither g nor r can be 
arbitrarily increased — both for fundamental and techni- 
cal reasons — , as we will show in Sec. II B. Note that no 
need not be equal to the mean bath occupation n, but 
may be decreased by laser pre-cooling to improve the 
protocol's performance. 

3. Entanglement verification 

To verify the successful creation of entanglement a red 
detuned laser pulse (A c = uj m ) is sent to the cavity 



where it resonantly drives the beam-splitter interaction, 
and hence generates a state swap between the mechanical 
and the optical mode. It is straightforward to show that 
choosing A c = w m leads to a different set of Langevin 
equations which can be obtained from (4) by dropping 
the Hermitian conjugation (f) on the right-hand side. 
By defining modified mode functions aJ n ( out ) = ctout(in) 
and corresponding light modes An(out) one obtains in- 
put/output expressions in analogy to (7) 

A' out = - c- Gt A[ n + iy/l - c- 2 Gr Bin , (10a) 

B ou t = e~ Gr B in - iy/l - e- 2 Gr A [ n . (10b) 

The pulsed state-swapping operation therefore also fea- 
tures an exponential scaling with Gr. For Gr — > oo the 
expressions above reduce to A' out = — iB m and B out = 
iA[ nl which shows that in this case the mechanical state — 
apart from a phase shift — is perfectly transferred to the 
optical mode. In the Schrodinger-picture this amounts 
to the transformation |<p) m | , 0)i — > \ip) m \tp)i, where tp and 
i\) constitute the initial state of the mechanics and the 
light pulse, respectively. The state-swap operation thus 
allows us to access mechanical quadratures by measur- 
ing quadratures of the light and therefore to reconstruct 
the state of the bipartite system via optical homodyne 
tomography. For this the protocol is operated in two 
steps: After the blue detuned pulse is reflected from 
the cavity, it is sent to a homodyne detection setup, 
where a quadrature X\(4>) — Aj out cos0 + Bj out sin0 (<j) 
being the local oscillator phase) is measured. The same 
procedure is subsequently carried out for a red detuned 
pulse, measuring X[{6) (see Fig. 1), which, for the case 
of a perfect state swap, yields the mechanical quadrature 
X[(6) = X m (0+§). Here the rotation by f is due to the 
phase shift from the swap operation. By repeating this 
process multiple times for different local-oscillator phases 
(4>i, 9j), the quantum state of the bipartite optomechan- 
ical system can be reconstructed. Having obtained the 
full quantum state, entanglement can be analyzed by var- 
ious means [53], e. g., by applying the EPR criterion from 
above. 



4- Optomechanical teleportation protocol 

As we have shown above, pulsed operation allows us 
to create EPR-type entanglement, which forms the cen- 
tral entanglement resource of many quantum information 
processing protocols [41]. An immediate extension this 
scheme is an optomechanical continuous variables quan- 
tum teleportation protocol. The main idea of quantum 
state teleportation in this context is to transfer an arbi- 
trary quantum state \ipin) of a traveling wave light pulse 
onto the mechanical resonator, without any direct in- 
teraction between the two systems, but by making use 
of optomechanical entanglement. The scheme works in 
full analogy to the CV teleportation protocol for pho- 
tons [54, 55]. Due to its pulsed nature it closely resem- 
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bles the scheme used in atomic ensembles [35, 56] and 
it was recently also suggested in the context of levitated 
microspheres [36] (see [57] for an exhaustive description 
of a similar system comprising a nuclear-spin ensemble 
entangled with light): A light pulse (A) is sent to the op- 
tomechanical cavity and is entangled with its mechanical 
mode (B) via the dynamics described above. Meanwhile 
a second pulse (V) is prepared in the state |"0i n ), which is 
to be teleported. This pulse then interferes with A on a 
beam splitter. In the output ports of the beam splitter, 
two homodyne detectors measure two joint quadratures 
pout _j_ Xv and X out + y i e idi n g outcomes m x and 

mp, respectively. This constitutes the analogue to the 
Bell measurement in the case of qubit teleportation and 
effectively projects previously unrelated systems A and 
V onto an EPR state [58]. Note that both the second 
pulse and the local oscillator for the homodyne measure- 
ments must be mode-matched to A after the interaction; 
i.e., they must possess the identical carrier frequency as 
well as the same exponential envelope. The protocol is 
concluded by displacing the mirror in position and mo- 
mentum by mx and mp according to the outcome of 
the Bell-measurement. This can be achieved by means 
of short light pulses, applying the methods described in 
[43, 45]. After the feedback the mirror is then described 
by [41] 

Xfin x^out i r>out , v" 
m — A m +M + A v. 

- A v + (> -y/c 2r -l) - Ph, (Ha) 

Pifin rjout i \rout i / > 
m — r m +A[ +r v , 

= P V + (e r - y/tfr -l) (i£ - Xl a ) , (1 lb) 

which shows that its final state corresponds to the in- 
put state plus quantum noise contributions. It is ob- 
vious from these expressions that the total noise added 
to both quadratures [second term in (11a) and (lib), re- 
spectively] is equal to the EPR variance. Again, for large 
squeezing r » 1 the noise terms are exponentially sup- 
pressed and in the limit r — > oo, where the resource state 
approaches the EPR state, we obtain perfect teleporta- 
tion fidelity, i.e., X^ 1 — X v and P^ 1 = P v . In particular 
this operator identity means, that all moments of X v , P v 
with respect to the input state \ipin) will be transferred 
to the mechanical oscillator, and hence its final state will 
be identically given by IV'in)- 

To verify the success of the teleportation one has to 
read out the mirror state after completing the feedback 
step. This can be achieved by applying tomography on 
the mechanical state as described in the previous section. 
The overlap of the reconstructed state p ut and a pure 
input state \ipi n ) then gives the teleportation fidelity F = 
(V'inl/foutlV'in)- For coherent input states the fidelity is 
given by 

F=(i+*f*y\ (12) 



In order to beat the optimal classical strategy for trans- 
mission of quantum states (i.e., the optimal measure- 
and-prepare scheme) , the achieved fidelity (averaged over 
all coherent states) must exceed F > 1/2 [59], which is 
equivalent to the condition for entanglement, Aepr < 2. 



B. Optimised protocol including imperfections 

1, Perturbations 

In the previous section we found that in the ideal sce- 
nario the amount of entanglement essentially depends 
only on the coupling strength (or equivalently on the in- 
put laser power) and the duration of the laser pulse and 
that it shows an encouraging scaling, growing exponen- 
tially with Gt. This in turn means that the minimal 
amount of squeezing needed to generate entanglement 
only grows logarithmically with the initial mechanical 
occupation hq. In this section we will develop a more 
realistic scenario including thermal noise effects and full 
system dynamics, both of which will decrease the cre- 
ated entanglement. We will show, however, that under 
conditions already available in state-of-the-art optome- 
chanical experiments, one can find an optimal working 
point such that the significance of these unwanted effects 
can sufficiently be suppressed. 

To extend the validity of the previous, simplified 
model, we now include the following additional dynam- 
ics: contributions from the beam-splitter Hamiltonian, 
higher order interactions beyond the adiabatic approxi- 
mation, and decoherence effects due to mechanical cou- 
pling to a heat bath. In the following we will investigate 
their effect on our protocol and determine the parameter 
regime featuring maximal entanglement. The technical 
details of how we include them in our calculations will be 
shown in Sec. III. 

Including the above-mentioned perturbations results in 
a final state which deviates from an EPR-entangled state. 
To minimise the extent of these deviations, the system 
parameters must obey the following conditions: 

1. k <C w m results in a sharply peaked cavity response 
and implies that the down-conversion dynamics is 
heavily enhanced with respect to the suppressed 
beam-splitter interaction. 

2. g < k inhibits multiple interactions of a single pho- 
ton with the mechanical mode before it leaves the 
cavity. This suppresses spurious correlations to the 
intracavity field. It also minimises pulse distortion 
and simplifies the protocol with regard to mode 
matching and detection. 

3. gr ^> 1 is needed in order to create sufficiently 
strong entanglement. This is due to the fact that 
the squeezing parameter r = (g/n) gr should be 
large, while g/n needs to be small. 



6 



4. njr <C 1, where n is the thermal occupation of the 
mechanical bath, assures coherent dynamics over 
the full duration of the protocol, which is an es- 
sential requirement for observing quantum effects. 
As the thermal occupation of the mechanical bath 
may be considerably large even at cryogenic tem- 
peratures, this poses (for fixed 7 and n) a very strict 
upper limit to the pulse duration r. 

Note however that not all of these inequalities have 
to be fulfilled equally strictly, but there rather exists an 
optimum which arises from balancing all contributions. 
It turns out that fulfilling (4) is critical for successful 
teleportation, whereas (l)-(3) only need to be weakly 
satisfied. Taking the above considerations into account, 
we find a sequence of parameter inequalities 

n-y < - < g < k < w m , (13) 

T 

which defines the optimal parameter regime. In Sec. II A 
we assumed the first two conditions to be well satisfied 
and we neglected the existence of mechanical decoher- 
ence. If we now take into account that the mechanical 
oscillator couples to a heat bath with an effective deco- 
herence rate 717, we find that increasing the pulse dura- 
tion to values larger than the mechanical coherence time 
will drastically decrease entanglement. This results in an 
upper bound for entanglement, as now both the interac- 
tion strength and the pulse duration, and therefore also 
the squeezing parameter r = (g/n)gT, are bounded from 
above. 

Dividing (13) by 7 and taking a look at the outermost 
condition n <C Q m , where Q m — w m /7 is the mechanical 
quality factor, we see that the ratio Q m /n defines the 
range which all the other parameters have to fit into. It 
is intuitively clear, that a high quality factor and a low 
bath occupation number, and consequently a low effec- 
tive mechanical decoherence rate, are favourable for the 
success of the protocol. Equivalently, we can rewrite the 
occupation number as n = fceTbath/Swrn and therefore 
find fceXbath/fi- ^ *9m ' w mj wnere now the Q 1 /-product 
(/ = w m /27r) has to be compared to the thermal fre- 
quency of the bath. Let us consider a numerical exam- 
ple: For a temperature Tbath ~ 100 mK the left-hand side 
gives /ceTbath/^ ~ 27r • 10 9 Hz. The Q ■ /-product con- 
sequently has to be several orders of magnitude larger 
to successfully create entanglement. As current optome- 
chanical systems feature a Q-/-product of 2ir- 10 11 Hz and 
above [60-63], this requirement seems feasible to meet. 
Note that in an experiment Tbath will often depend on 
the input laser power, as scattered light can heat up the 
cryogenic environment. Hence, for a given bath occupa- 
tion the coupling strength may be limited for technical 
reasons. 

In order to find the optimal working point, it is con- 
venient to introduce the following dimensionless parame- 
ters: the sideband-resolution parameter rj, the adiabatic- 
ity parameter £, and the ratio of pulse length to mechan- 



ical coherence time e. They are given by 

V = K/u) m , £ = g/n, e = 77-. 

From (13) it follows that rj < 1, £ -C 1 and e < l/n. 
Each of those small parameters can be used to realise a 
perturbative expansion of the additional dynamics listed 
above. The perturbative solutions can then be used to 
calculate the EPR variance and optimise the resulting 
expressions. 

2. Optimization 

As illustrated above, we expect — for fixed values of n 
and Q m — to find optimal values for the remaining param- 
eters e, £ and rj. The maximal possible entanglement will 
ultimately be set by n and Q m , which will also constitute 
hard boundaries in typical experiments. 

Figure 2 shows results of this optimization for different 
Q-factors, plotted against the thermal occupation num- 
ber of the mechanical bath. Figure (a) shows the min- 
imal value of Aepr for a given Q m , and (b)-(d) show 
the corresponding optimal values for rj op t 7 £ opt and e op t. 
As expected the noise contribution from the mechanical 
bath is found to be the most critical. As we show in 
IIIB the EPR variance including thermal noise can be 
expressed by 

A EPR - [A(X^ + Pr')] 2 

+ [A(P° ut + ^i° ut )] 2 + (2n + l)e, (14) 

where X° ut and P° ut here denote the solutions for 7 = 0. 
The dashed curves in Fig. 2 (a) illustrate the noise contri- 
bution (2h+ l)e op t of the thermal bath. As shown above, 
this quantity is added to the unperturbed EPR variance, 
and thus the system can only exhibit strong entanglement 
if its value is far below two. Note that working at the op- 
timal point keeps the fraction of thermal noise in Aepr 
approximately constant over a wide range of h. This is 
shown in Fig. 2 (e), where we defined e opt = e op t/AEPR. 

As we have seen in the previous section, the EPR vari- 
ance depends on the occupation number of the oscillator 
at the initial time t = 0. Due to this, the entanglement 
can be drastically increased by pre-cooling the mechanics 
by means of laser cooling before starting the actual pro- 
tocol. Fig. 2 shows that it is thus possible to create an 
entangled state even for a fairly large bath occupation. 
This works due to short pulse durations, during which 
the mechanical decoherence is small. 

Taking a look at figure (c) we note that the sideband- 
resolution shows rather large optimal values near unity, 
especially for increasing occupation numbers. This indi- 
cates that the beam-splitter dynamics only weakly dis- 
turbs the entangling interaction. 

Table I gives a list of the optimized key experimental 
parameters for two existing optomechanical structures. 
The values in the first row correspond to the black dots 
shown in Fig. 2. 
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(a) 

Tbath [mK] at 3.8MHz 

0.0 45.7 91.3 137.0 182.0 228.0 274.0 
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(b) (c) 




500 1000 1500 500 _ 1000 1™° 



n n 

FIG. 2. (Color online) Optimised parameters for Q m = 10 7 
and no = 50 [red (dark) line], and Q m = 10 and no = [blue 
(light) line], where no is the initial mechanical occupation and 
n is the mean bath occupation. This corresponds to the two 
cases of large Q m with moderate pre-laser-cooling, and lower 
Qm with pre-cooling into the ground state. Clearly, entangle- 
ment creation is possible in both cases, (a) Minimal Aepr 
as a function of n. The upper axis gives the corresponding 
bath temperature for a oscillator with a resonance frequency 
of 3.8 MHz. To each point corresponds a triple e op t , r) opt , £ opt 
[(b)-(d)] for which the minimal value is realised. The dotted, 
black line shows the upper bound up to which entanglement 
is present. The black dots marks the values for n = 1100 
(T = 200 mK), which coincide with the values given in the 
first row of table I. The dashed lines show the respective ther- 
mal noise contribution (2n+l)e opt to Aepr [see (14)]. (b)-(d) 
Values of e, rj, £ which optimise Aepr for a given n. (e) Shows 
the relative amount of noise induced by the coupling to the 
mechanical environment (e' opt = e op t/AEPR). 



Over and above the fundamental imperfections, techni- 
cal losses, such as mode- mismatch and detector inefficien- 
cies, additionally decrease entanglement. They can all 
collectively be described as (passive) beam-splitter losses 
adding vacuum noise to the optical signal. They will 
however, never completely break entanglement, as long 



as the overall loss is smaller than unity. Noise contribu- 
tions of this type can easily be accounted for by adding 
appropriate noise terms to (7). 

Finally let us compare the amount of entanglement 
created in the pulsed and continuous- wave schemes [14] 
in terms of the logarithmic negativity £jv [53] for the 
parameters used in the first row of table I (Q m = 10 5 , n = 
1100, no = 0). In the case of continuous driving one finds 
the maximal negativity close to the instability region and 
for a detuning of around A c sa w m , yielding E_\f sa 0.4. 
For the pulsed protocol the optimisation yields a much 
larger value of Ejs/ ss 1.2. Note that Ejs/ is a logarithmic 
quantity. 



III. DETAILED MODEL 
A. Linearizing the dynamics 

As was shown in [65], the radiation pressure interac- 
tion is inherently non-linear. In current micromechanical 
systems however, the single-photon coupling g$ is very 
weak (the best values up to date are on the order of 
go/ft 0.001 [63]), and has to be enhanced by means 
of a strong optical pump field. It is well known that 
in the case of a strong continuous-wave light field, the 
steady-state dynamics of the system is approximately lin- 
ear. We will show in the following, that this also holds in 
a (long) pulsed scheme. The linearization process follows 
the same general idea as in the steady-state regime, it is, 
however, slightly more involved due to the explicit time 
dependence of the Hamiltonian. 

We consider a laser pulse with a fixed number of pho- 
tons N p h and an envelope function e(t), which is nor- 
malised in the sense that JJidi |e(i)| 2 = 1. Its head and 
tail are assumed to be smooth and its amplitude should 
be constant e(t) « 1/\/t for the most part of r. The full 
Hamiltonian for the system, including the laser driving 
term and the non-linear radiation pressure interaction 
[65] is then given by 

H(t) = w m o4a m + A ala c + 

g ata c (a m + a f m ) + \E{t) (a c - aj) , (15) 

where Ao = w Cj o — is the detuning for the case of 
a cavity with fixed length and E(t) = yj2nN- p \ l e{t) is 
the driving strength. In the Appendix we show that we 
can eliminate the driving term by going into a (time- 
dependent) displaced picture. The transformed Hamil- 
tonian then takes the form 

H(t) = u m a\ a a m + \{t)a\a c + g a a\a c {a m + 4j 

+ 00 ("* OK + °<(t)ai) (a m + a L)> (16) 

where the effective detuning A c and the mean cavity field 
a now depend explicitly on time. The resulting expres- 
sions (see the Appendix) are essentially the same as found 
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TABLE I. Specific optimal values for an ALGai-^ As structure (3.8 MHz) [64] and an Si optomechanical crystal (3.7GHz) [10]. 
The whole set is fully determined by e op t, ^o P t, £o P t- The value for the mean power P = hcuiNph/r is obtained using (2). (The 
x in the formula above is a number between and 1 to indicate a ternary alloy between GaAs and AlAs.) 



CJ m /27T 


Qm 


Tbath 


n 


n 


g /2n 


K opt /27r 


T~opt 


Popt 


ffopt/27T 


Aepr 


3.8 MHz 


10 b 


200 mK 


1100 


0.0 


4.8 Hz 


3.2 MHz 


2.5 /is 


30 mW 


0.97 MHz 


0.7 


3.7 GHz 


10 5 


200 mK 


0.7 


0.7 


910.0 kHz 


0.26 GHz 


0.41 ^s 


6/xW 


0.032 GHz 


0.1 


3.7 GHz 


10 5 


IK 


3.7 


3.7 


910.0 kHz 


0.31GHz 


0.30 /is 


8/xW 


0.040 GHz 


0.5 



in the steady-state case (see for example [66]). The non- 
linear term in (16) can be neglected whenever |ev| 3> 1, 
which is true for sufficiently strong driving \E\ and will 
be the case for the greatest part of the pulse duration. 
By assuming a(t) to be real and by introducing the ef- 
fective optomechanical coupling constant g(t) = go \a(t)\, 
the procedure leaves us with a linear Hamiltonian in the 
form of (1). 

Note that for the case of A c w oj m the relative fre- 
quency shift induced by radiation pressure will be of the 
order of O(g/oj m ) 2 , and therefore small. Consequently 
[together with the assumption that e(t) s» l/y^r] we W1 U 
in the following drop the explicit time dependence of the 
effective detuning and the effective coupling strength. 



B. Solving the full system 

The full Langevin equations (3) resulting from the lin- 
earised Hamiltonian (including the beam-splitter inter- 
action and mechanical decoherence) can be rewritten in 
the compact form 



dt 



R(t) = (S- D)R(t) - V2DR in (t), 



(17) 



where R — (a m , d^, a c , aj) and correspondingly i?i n de- 
notes the input noise. S and D are matrices compris- 
ing the respective coefficients. To solve this set of equa- 
tions we apply the Laplace transformation, introducing 
R(s) = C[R](s), with £[/](s) 
for R we obtain 



R(s) = M(s) (r(0) 



and thus 



R(t) = M(t)R(0) - (m * 



(18) 



(19) 



where M(s) = (si - S + D)' 1 , M(t) = £" x [M](r) 
and * denotes the convolution integral (/ * g)(t) = 

J ds g(s)f(t — s). In the case of a bipartite system it is 
possible to find an exact expression for M for arbitrary 
parameters. The obtained solution, however, is very te- 
dious and will not be presented here. We proceed as 
follows: We separate the mechanical decoherence in a 
perturbative approach (i.e., we expand M in powers of 
e), while the other dynamics will be treated exactly. This 



allows us to find input /output relations corresponding to 
(7) and to calculate the EPR variance for the full system. 

We established in Sec. II B 2 that for the protocol to 
work we require the effective mechanical decoherence 
time to be much larger than the duration of the light 
pulse, i.e., en <C 1. We emphasise that e is the smallest 
of all parameters, and the coherent evolution will only 
be negligibly perturbed by the coupling to the mechani- 
cal bath. We will therefore only keep terms 0(en) while 
neglecting 0(e). Based on this premise we simplify (19) 
twofold: Firstly we drop the mechanical damping from 
the first term, as it gives corrections on the order of O(e) 
only; thus M w Af| 7=0 - This amounts to dropping the 
term -|p m in (3b). Secondly, we approximate the me- 
chanical noise contribution (second term) by only keep- 
ing the free, harmonic evolution, while neglecting their 
coupling to the optical mode. This coupling is due to a 
second-order process, and is therefore suppressed by an 
additional factor of £. Note that by doing so we overes- 
timate the effect of mechanical noise, as it contributes to 
the creation of optomechanical correlations when subject 
to the coherent dynamics. The complete noise term en- 
tering in the evolution of the mechanical variables then 
takes the form 



iy/j / ds e 
/o 



7(r 



(F s + iF c ) , (20) 



where we introduced (co)sine components F( c j s of the 
Brownian force. We can therefore write a m (r) w 



f^dt e~ st f(t). Solving «mW! 



7=0 



iF c ), and consequently 



^>out ^ pout 



7=0 
7=0 



+ VeF r ., 



(21a) 
(21b) 



while we neglect the mechanical noise contribution to the 
optical mode, i.e., A out ~ v4 out | 7= o- Note that from 
the commutation relation of the Brownian noise term, 
[f{t + a), f{t)] = ±6'(s) [50], it follows that [F„F e ] = 
i + 0(l/w m T) and [F^F,] = 0(l/w m r). The perturbed 
variables (21) therefore approximately obey canonical 
commutation relations [X^ ut , P™'] rj i(l + e). Using (21) 
together with the correlation functions (FiFi) = n + \ 
leads to (14). As we have separated the mechanical noise 
terms from the other dynamics, we will always assume 
that M w M| 7=0 an d drop the 7 dependence for the rest 
of this section. 
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We now use (19) together with the definitions of A ou t 
(6b) and B out to obtain input /output equations similar 
to (7), but for the full system dynamics. The resulting 
expressions are of the form 



Cl B in + c 2 Bj n + c 3 a c (0) + c 4 aj(0) 



+ c 5 / dsa in! i(s)a in (s) + c 6 / ds a* n2 (s)al n (s), 
Jo Jo 

with a similar expression for A out - The coefficients Ci as 
well as the light modes a- ln ^ are determined by the sys- 
tem dynamics [given by M(t)\ and the light mode a out 
selected from the output field. Note that these expres- 
sions are valid for 7 = only and thus have to be used in 
conjunction with (21) to account for mechanical noise. 

Following the treatment in Sec. II A 2 one easily finds 
the corresponding EPR variance, which now includes 
noise terms from the initial intracavity field and the extra 
light modes (both assumed to be in vacuum). The lat- 
ter contributions are given by the overlap of the different 
light modes /Jdt a in>i (t) a* lni (t). 

The resulting expression is an involved function of e, r\ 
and £ and is not presented here. Numerical minimization 
with respect to those three variables for fixed values of n 
and Q m yields the results presented in Fig. 2. 



IV. CONCLUSIONS 

We have developed a scheme to create and — due to 
its pump-probe operation — unambiguously verify EPR 
entanglement in optomechanical systems. Additionally 
its application as an entanglement resource in quantum 
teleportation was discussed. Finally, by optimizing the 
experimental parameters we showed that the suggested 
protocol is feasible with state-of-the-art optomechanical 
devices. 
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Appendix A: Transformed Hamiltonian 

Starting from the Hamiltonian (15) we write down the 
standard quantum-optical master equation for a damped 
cavity mode 



p = -i [H, p]+K (2a c pal - a c a\p - pa, 



(Al) 



while we neglected the mechanical decoherence terms as 
we are only interested in times far within the coherence 
time of the oscillator. In order to eliminate the driving 
field E(t) we go into a displaced picture 



p = D c (a)D m (p)pDl(a)Dl(l3), 



(A2) 



with displacement operators Di(a) — exp(aa| + a*ai). 
The time-dependent, complex amplitudes a = a(t) and 
(3 — f3(t) give the mean displacements due to the laser 
drive and will be determined in the following. 

The transformed master equation can again be written 
in the form of (Al) by substituting p — > p and H — » H . 
The Hamiltonian H is then given by 

H = u m a] a a m + [A + g (P + j3*)]a\a c 

+ g a (a*a c + aa\)(a m + ajj + g ala c (a m + ajj 

+ |[id+ (i« + A )a + g {/3 + /3*)a - \E\a\ +h.c.| 



i/3 + u m j3 + go\a\ 



(A3) 



The first two lines constitute the new Hamiltonian of the 
system, and the last two lines describe the mean (classi- 
cal) cavity and mirror amplitude, respectively. We can 
make these terms disappear by choosing a and ft such 
that they fulfill the following set of coupled, non-linear 
differential equations 

a = {i[A +g Q (p + p*)]-K}a + E, (A4a) 
P = iu m /3 + ig a \a\ 2 . (A4b) 

We seek solutions to these equations for initial conditions 
a(0) = (3(0) = in terms of the driving field E(t). Due 
to their non-linear nature, no exact closed-form solution 
will exist in general and we will therefore look for ap- 
proximate solutions under the assumptions we made in 
Sec. II B. We formally integrate equation (A4b) to find 



P(t)=igo / ds e^ s \a(t- 



(A5) 



Under the assumption that k -C w m and given that E(t) 
varies sufficiently slowly, we also expect |a(t)| to be a 
slowly varying function on the time scale of l/cJ m . We 
will check this for consistency at the end of this section. 
For this case we can use the adiabatic solution 

9o 



Ul n 



(A6) 



Plugging this into (A4a) and introducing the effective 
detuning 



Ac(t) = Ao-floG0(*) +£*(*)) 
25o 2 



A 



„ \a(t)\' 



we find the solution 



1 _ p -(iA c +re)t 



E(t) 

lA r + K ' 



(A7) 



(A8) 
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where 5 is a correction, which is small if E(t) varies slowly 
on a timescale of 1 /k. More precisely, one can show that a 

crude upper bound is given by \5(t)\ < sup se(o t) ~j§^[, 
which must be much smaller than unity. Also, as we as- 
sume that kt > 1, we neglect the term e~ K *, as this only 
contributes at the very beginning of the pulse. The ap- 
proximations made in (A8) amount to assuming that the 
slope of the pulse is small enough (with respect to k) that 
it does not experience distortion due to the finite cavity 
linewidth. Throughout this derivation we have assumed 
that A c (t) w A c is approximately constant in time, which 



is well fulfilled for the parameter regime O(g/oj m ) 2 <C 1 
that we are concerned with. This can be made exact if 
we assume the laser detuning to be locked with respect 
to the effective cavity resonance frequency. Having ob- 
tained a solution for a(t), we can go back to test the 
self-consistency of our derivation of (A6), where we re- 

quired U <C w m . In the case g <C K -C w m we find 

t ha t ^l^^ml^ 

w m , giving an additional condi- 
tion on the pulse shape. Taking all these considerations 
into account we arrive at the linearised Hamiltonian (16). 
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